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ABSTRACT 

The  Ziv-Zakai  bound  (ZZB)  provides  a  general  mean-square  error 
analytical  baseline  to  evaluate  time  delay  estimation  (TDE)  tech¬ 
niques  for  a  wide  range  of  time-bandwidth  products  and  signal-to- 
noise  ratios,  but  generally  can  only  be  numerically  evaluated.  The 
Weiss-Weinstein  bound  (WWB)  further  improves  characterization 
of  the  attainable  system  performance,  for  narrowband  and  wide¬ 
band  signals  with  small  to  moderate  fractional  bandwidth.  Sim¬ 
ilar  to  the  WWB,  here  we  find  a  simplified  closed-form  ZZB  for 
TDE  with  ultra-wideband  (UWB)  signals.  The  resulting  simplified 
bound  is  found  over  disjoint  segments,  separated  by  thresholds  that 
characterize  different  regions  of  ambiguity.  The  closed-form  sim¬ 
plified  bound  can  be  analytically  studied,  and  approaches  both  the 
ZZB  and  the  TDE  performance  of  a  maximum  likelihood  estima¬ 
tor. 

INDEX  TERMS 

Ultra-wideband,  time  delay  estimation,  Ziv-Zakai  bound,  Weiss- 
Weinstein  bound. 

1.  INTRODUCTION 

Time  delay  estimation  (TDE)  is  fundamental  to  many  applications, 
such  as  array  processing,  localization,  and  tracking.  Time  delay 
estimator  performance  is  often  characterized  via  mean  square  er¬ 
ror  (MSE),  that  can  then  be  compared  to  analytical  MSE  lower 
bounds.  The  Cramer-Rao  bound  (CRB)  predicts  maximum  likeli¬ 
hood  estimator  (MLE)  performance  for  sufficient  time  bandwidth 
(TB)  product  and  signal  to  noise  ratio  (SNR)  [1].  However,  the 
CRB  is  a  local  bound  that  fails  to  characterize  performance  when 
the  ambiguity-free  condition  is  violated,  for  example,  when  a  cross¬ 
correlation  TDE  incorrectly  selects  neighboring  correlation  peaks. 
Tighter  bounds  for  TDE  have  been  developed,  including  the  Chapman- 
Robbins  (Barankin)  bound  [2],  and  the  Ziv-Zakai  bound  (ZZB) 

[3].  Though  much  tighter  than  the  corresponding  CRB,  they  are 
not  easily  manipulated  into  simple  closed  form  expressions,  and 
consequently  often  require  numerical  evaluation.  Eurther  improve¬ 
ments  have  been  made  by  Bellini  &  Tartara  [4],  and  Weiss  &  Wein¬ 
stein  [5],  [6].  Extension  to  the  vector  case  has  also  been  developed 
[7].  The  improved  ZZBs  in  [5],  [6]  are  applicable  to  either  nar¬ 
rowband  or  wideband  waveforms  because  of  explicit  assumptions 
made  throughout  the  derivations.  This  facilitates  closed  form  ex¬ 
pressions  in  some  cases,  or  close  approximations  that  have  simple 
forms.  TDE  for  frequency  hopping  signals  is  considered  in  [8], 
and  an  overview  of  TDE  bounds  is  given  in  [9]. 

Ultra-wideband  (UWB)  signals  have  long  been  employed  in 
radar,  and  more  recently  for  communications  and  geolocation,  e.g.. 
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in  the  3.1GHz  to  10.6GHz  band  in  the  US  [10,  11].  Timing  ac¬ 
quisition,  synchronization,  and  delay  estimation  are  important  and 
challenging  aspects  of  these  systems,  with  sensitivity  to  jitter  and 
signal  mismatch  [12,  13].  We  consider  TDE  in  the  UWB  regime, 
assuming  large  fractional  bandwidth,  and  simplify  the  ZZB  for  this 
case.  Similar  to  the  developments  in  [5,  6]  for  the  narrowband  and 
wideband  cases,  we  develop  the  ZZB  for  the  UWB  case,  and  find  a 
piece-wise  approximation  in  simple  closed  form  that  is  amenable 
to  analysis  in  applications.  Eor  the  UWB  case,  the  correlation  ex¬ 
hibits  a  significantly  different  pattern  than  for  signals  with  smaller 
fractional  bandwidth,  such  that  there  is  no  clearly  defined  phase  or 
envelope  ambiguity  corresponding  to  different  correlation  peaks. 
Consequently,  we  develop  a  new  interpretation  of  ambiguity,  that 
depends  on  the  lower  and  upper  signal  frequencies.  Similarly,  we 
define  new  threshold  points.  The  new  ambiguities  arise  in  seg¬ 
ments  dominated  by  the  envelope,  phase  and  their  mixtures,  that 
in  turn  depend  on  TB  product  and  SNR.  Einally,  a  simulation  study 
shows  that  our  simplified  ZZB  closely  predicts  MLE  performance, 
and  approaches  the  true  ZZB  (obtained  numerically)  in  the  differ¬ 
ent  regions. 

2.  REVIEW  OE  ZZB  AND  WWB 

Consider  estimation  of  delay  r  between  two  sensors.  Each  wave¬ 
form  is  cormpted  by  additive  noise  so  that 

ri(f)  =  s{t)+ni{t), 

r2{t)  =  s{t-T)  +  n2{t)  -Tl2<t<Tl2,  (1) 

where  T  is  the  observation  time,  and  s{t),  ni{t)  and  n2{t)  are 
the  sample  functions  of  the  signal  and  Gaussian  noise  with  spec¬ 
tral  densities  S{u!),  Ni{uj),  and  ^2(0;),  respectively.  The  ZZB 
development  begins  with  a  binary  decision  problem  based  on  an 
arbitrary  estimate  r  of  r. 

Decide  Ho  ■.  t  =  a  if|r  —  a|<|r  —  a  —  0|, 

Decide  Hi  ■.  t  =  a  -\-  6  if|r  —  a|>|r  —  a  —  0|.  (2) 

Without  loss  of  generality,  the  two  hypothesized  delays  are  as¬ 
sumed  to  be  equally  likely  to  occur,  so  the  probability  of  error  is 
P{t  — a  >  0/2\t  =  a\/2  +  P{T  —  a  —  6  <  —9/2\t  =  a  +  9}/2. 
Denote  the  error  in  estimating  t  hyr  os  t  =  T  —  t. 

We  are  interested  in  the  mean  square  error  =  E{te^)  with 
T  uniformly  distributed  within  the  interval  [—D/2,  D/2],  and  9 
within  [0,  Dj.  The  basic  ZZB  is  given  by  [3] 

—  1  rD/2-e 

>  -—  /  9  d9  Pe{a,a  +  9)da,  (3) 

D  Jo  J -D/2 
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where  Pe(a,  a  +  6^)  is  the  minimum  attainable  probability  of  error 
associated  with  the  likelihood  ratio  test  between  the  two  hypothe¬ 
sized  delays.  It  immediately  follows  that 

Pe{a,a+B)  <  >  S/2\t  =  a)-|-P(e  <  —0/2\t  =  a-|-6»)}. 

(4) 

Bellini  and  Tartara  further  noticed  that  ^  ci+d)da 

is  a  nonincreasing  function  of  0  [4].  Therefore,  by  using  a  nonin¬ 
creasing  function  V[-]  that  fills  the  valleys  of  the  function  inside 
the  bracket,  an  improved  ZZB  is  given  by 


ev 


l<D/2-0 


Ps{a,  a  +  0)  da 


-D/2 


de.  (5) 


If  Pe{a,a  -f  0)  =  P<i{S)  is  independent  of  a,  the  improved 
ZZB  may  be  simplified  to  the  WWB  [5] 

ev[{D-e)p^{e)]de.  (6) 

Typically,  the  observation  time  T  is  large  compared  to  the  cor¬ 
relation  time  2-k  jW  where  II^  is  the  bandwidth.  Then,  Pe(0)  is 
closely  approximated  by 


P^[e)  «  e“(«)+''(»)<I>(y^2fo(0)),  (7) 

where  $(*)  =  ^  df,. 


P  OO 

a{9)  = - /  ln[l  -f  SNR(t<;)  sin^(t<;0/2)]  duj, 

2^  Jo 


(8) 


m 


T  SNR(a;)sin2(w6»/2) 

^  H-SNR(w)sin2(w6»/2) 


(9) 


=  [S{o2)/N^{lo)][S{lo)/N,{lo)] 

^  ’  1  +  S{lo)/Ni{uj)  +  S{uj)/N2{uj)' 


(10) 


In  the  derivation  of  (7),  large  WT  is  assumed.  This  condition  is 
satisfied  in  most  practical  scenarios.  Although  it  is  difficult  to  sim¬ 
plify  (8)  and  (9)  in  general,  [5]  and  [6]  considered  narrowband 
and  wideband  waveforms,  and  derived  closed  form  expressions 
for  different  disjoint  segments  separated  by  thresholds.  The  cor¬ 
responding  simplified  ZZBs  are  very  valuable  for  comprehensive 
analytical  study  of  bound  performance  under  different  system  op¬ 
erating  conditions  when  the  fractional  bandwidth  of  the  waveform 
is  small  or  moderately  large  (the  “narrowband”  and  “wideband” 
cases).  In  the  following,  we  develop  a  simplified  and  analytically 
tractable  ZZB  for  the  UWB  case  with  large  fractional  bandwidth. 


3.  SIMPLIFIED  ZZB  FOR  UWB  SIGNALS 

We  consider  source  signals  and  noises  that  are  spectrally  flat  over 
the  frequency  band  [wo  —  W/2,tjJo  -f  IU/2]  with  spectral  densi¬ 
ties  S  and  N  respectively,  where  uq  denotes  the  center  frequency. 
Cases  for  signals  transmitted  through  frequency  selective  multi- 
path  channels  will  be  studied  in  the  future.  Then,  the  SNR  defined 
in  (10)  is  constant  with  respect  to  W.  We  define  loh  =  uo  -f  IU/2, 
UJL  =  Wo  -  W/2,  rj  =  WT  ■  SNR,  and  A  =  SNRsin^(w6»/2) 
(this  last  term  appeared  in  (8)  and  (9)). 


The  most  common  estimator  of  r  is  to  find  the  peak  of  the 
cross-correlation  of  ri(f)  and  r2(t).  In  the  absence  of  noise,  this 
reduces  to  the  signal  auto-correlation 

Rss(j)  =  —  sinc(wirr)  —  —  sinc(wLr) 

TT  TT 

w  ^  ^ 

=  — cos(tJor)sinc(Wr/2).  (11) 

TT 

Although  the  peak  occurs  at  the  true  delay  r,  there  are  many  other 
local  maxima  t  =  t  -\-  9  that  appear  quasi-periodically.  In  the 
narrowband  and  wideband  cases,  typically  wo  W .  Then,  ac¬ 
cording  to  the  second  equality  of  (11),  phase  and  envelope  ambi¬ 
guities  occur  at  points  approximately  uniformly  distributed  with 
period  2-k  jwo  and  2k  jW ,  respectively  [6].  However,  as  W  grows 
closer  to  ujo,  cosjcuor)  cannot  be  perceived  as  a  carrier  signal  and 
sinc(IUr/2)  as  the  envelope.  However,  the  first  equality  of  (11) 
shows  that  sinc(wj^r)  changes  more  rapidly  than  sincjtUir)  when 
wh  wl-  This  fact  permits  us  to  similarly  introduce  the  con¬ 
cept  of  phase  and  envelope  ambiguities  for  the  UWB  case.  The 
phase  ambiguity  points  are  caused  by  u> h sine {u>hk)/k  with  pe¬ 
riod  2k/ujh,  and  the  envelope  ambiguity  points  are  caused  by 
t<;Lsinc(aJi,r)/7r  with  period  2k /wl-  Fig.  1  demonstrates  this 
transition  behavior  of  Ras^Jf)  with  wo  =  201U,  wo  =  5IU  and 
uio  =  W  respectively. 


Fig.  1.  Rasijr)  for  narrowband,  wideband,  and  UWB  signals. 

In  order  to  obtain  an  analytically  tractable  expression  for  (6) 
suitable  for  UWB  signals,  we  will  partition  the  integration  interval 
[0,79]  into  sub-intervals  separated  by  si,  S2  and  S3,  defined  be¬ 
low.  These  points  will  designate  the  effects  of  individual  phase  or 
envelope  ambiguity,  or  their  mixture.  Then  the  valley-filling  func¬ 
tion  12 [(79  —  9)Pa{9)\  will  be  closely  bounded.  Subsequently,  an 
analytical  expression  for  the  integral  of  each  sub-interval  will  be 
derived. 

3.1.  Simplification  and  Segmentation  of  V[(79  —  9)Pa{9)] 

As  previously  explained,  the  more  critical  ambiguity  results  from 
the  phase  of  the  autocorrelation  function  because  of  its  oscilla¬ 
tory  nature  with  period  2k  juiH,  and  the  secondary  ambiguity  is 
caused  by  the  oscillatory  nature  of  the  envelope  of  the  autocorre¬ 
lation  function  with  period  2k  juiL-  Therefore,  for  9  that  is  not 
very  close  to  0,  a  simplified  bound  on  V[(79  —  9)Pe{9)]  can  be 
developed  from  the  local  maxima  points  (i.e.,  both  the  phase  and 
the  envelope  ambiguity  points)  of  Pe{9),  or  equivalently  the  local 
maxima  (ambiguity  points)  of  the  autocorrelation  function.  When 
0  <  9  <  si,  no  ambiguity  needs  to  be  considered  in  the  au¬ 
tocorrelation  function,  and  this  region  is  noise-dominated.  When 
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Si  <  6^  <  S2,  the  phase  ambiguity  dominates.  Whens2  <  0  <  S3, 
both  the  phase  and  envelope  peaks  of  the  autocorrelation  function 
contribute  to  the  ambiguity.  Finally,  in  the  S3  <  9  <  D  region, 
the  envelope  ambiguity  has  the  dominant  effect.  In  the  following 
we  present  the  key  results,  and  omit  the  derivations  due  to  space 
limitations. 

For  the  noise  dominated  region  (0  <  0  <  si),  the  valley¬ 
filling  function  V{{D  —  6)Pe{6)]  is  bounded  by 


V[{D-e)P49)]  >  (12) 

where 

ci{0)  =  9^-^{12ujI  +  W^), 

di{9)  =  04.  (13) 


In  the  phase  ambiguity  dominated  region  (si  <  9  <  S2), 
the  peaks  of  the  autocorrelation  function  are  nearly  periodically 
distributed,  and  located  approximately  at  9„  =  21111/ ujh{  n  = 
1,2,...).  At  each  0n-i  <  9  <  9„,  the  valley-filling  function  can 
be  closely  bounded  by 

V[(£»-0)Pe(6')]  >  (P- — 

luh 

■  ^[\/  C2{9  +  2n/ujH)],  (14) 


where 


C2{9) 

d2{9) 


=  9- 


LOoV 


27rVl  +  2SNR 
j  (0)  I  tlSNR^  6 

.(4  +  ^^oVV^a.oV"  +  ^).  (15) 


In  the  region  (S2  <  9  <  S3),  where  both  phase  and  envelope 
ambiguities  are  dominant,  the  ambiguity  points  are  the  autocorre¬ 
lation  peaks  not  only  caused  by  the  highly  oscillatory  nature  of  the 
phase  but  also  the  nearby  envelope.  Suppose  9o,  9i,  02,  ...,  are 
the  local  maxima  of  Pe{9),  which  are  phase  ambiguity  points  and 
are  in  the  adjacent  area  of  the  envelope  ambiguity  points  as  well. 
Then,  2n/uJH  <  9n  —  9n-i  <  2it/u>l-  The  valley-filling  function 
is  strictly  bounded  as 


V{{D  -  9)P.{9)]  >{D-  —  -9)-  (16) 


with 


C3 


2ujhT  ■  SNR 
ttVI  +  2SNR 


,d3  =  ^(^turrSNR"  +  A,^^SNR®). 

(17) 


Here,  both  C3  and  da  are  independent  of  0. 

In  the  envelope  ambiguity  dominated  region  (S3  <  0  <  D), 
01,  02,  ... ,  are  the  local  maxima  associated  with  the  envelope  am¬ 
biguity  with  9n  =  2rm /lol-  The  valley-filling  function  can  be 
closely  bounded  as 


r  (P-27r/cui, -0)e-‘*"-l-(VSI) 
V[{D-9)P,{9)]>  {  {9<D-2iiIujl)  (18) 

i  0  {9>D-  2n/u;L), 


where  C4  =  2b  and  d4  =  —  (a  -F  b),  and 


TW  ,  1  -F  v/l  +  SNR  , 

a  = - In - - - ,b^ 

IT  2 


Both  C4  and  di  are  independent  of  0. 


TW  yi  +  SNR  -  1 
27r  -F  SNR 

(19) 


3.2.  Simplified  ZZB 


Now  we  collect  the  above  results  to  formulate  a  simplified  ZZB, 
and  consider  some  limiting  cases.  Partitioning  the  integration  in¬ 
terval  [0,  D]  in  (6)  into  segments,  and  substituting  (12),  (14),  (16) 
and  (18)  into  (6),  the  simplified  ZZB  for  UWB  signals  is  given  by 


> 


-F 


D 

1 

D 


9{D  -  0)e-'^il®l4-(i/ci(0))  d9 


2it 

9{D-  —  -0)e 


4?(\/c2(d  -F  27r/wrf))  d0 


1 

r"  .  27r 

-F 

'd  ^ 

/  0{D - 

Is, 

0)  •  e 

1 

27r 

-F 

/ 

D 

Is3 

Ti>L 

(20) 


In  this  expression,  si,  S2  and  S3  are  found  as  follows. 

Since  in  the  range  [0,  si),  noise  dominates  without  any  ambi¬ 
guity,  we  have 

si  =  2it/ujh-  (21) 

At  0  =  S2,  the  second  term  in  (20)  is  approximately  equal  to  the 
third  term,  therefore 


C2(S2  -F  2it/ujh)  =  C3. 


(22) 


At  0  =  S3,  the  third  term  in  (20)  is  approximately  equal  to  the 
fourth  term.  Then  S3  can  be  found  from 


TSNR 


[IH-F 


iTyj2{l  -F  2SNR)  2(s3  -  2it/ujl) 


a. 


(23) 


With  these  values  for  si,  S2  and  S3,  the  lower  bound  (20)  can 
be  analytically  studied  for  different  system  parameters.  Next  we 
consider  several  limiting  cases. 

When  ri  ^  0,  the  fourth  term  in  (20)  dominates,  and  the  lower 
bound  is  close  to 


>  1 
-  D 


cD—2-rr/bJi, 


/  0{D  -  27r/LJL  -  0)d0  • 

lo 

^{\/v/‘2tt)  ~  (24) 


When  1  7]  8v/2(H7  -F  2u>o)/u!o,  i-e.,  when  S2  <  S3,  the 

third  term  in  (20)  dominates,  and  the  lower  bound  is  very  close  to 


> 


si 


^  .  g-da 


^  ^2, 


1  27r 

2  W  2y^rjW' 

4  WL  2y^Wr] 


2!Z) 

2  ^ 


(25) 


If  0  ^  8V2(W  +  2tuo)/^<.^o,  the  first  term  in  the  summation 
becomes  the  dominant  one,  and  the  lower  bound  approaches 


> 


1 

c  J 
1 

4c 


n2-rr,/c/uJH 


9e 


-eV(P/d) 


<f>(0)  d9 


TliLol  +  1HV12)  ’ 


(26) 
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where  c  =  ci(0)/0^  =  77(12010  +  VK^)/487r. 

The  transition  from  (24)  to  (25)  starts  at  the  3dB  lower  point 
rj  —  a,  which  satisfies  4>(y^a/27r)  =  i  and  completes  at  the 
point  where  the  fourth  term  in  (20)  is  equal  to  the  third  term,  i.e., 
77  =  /3  when 


[2]  E.  W.  Barankin,  “Locally  best  unbiased  estimates,”  Awn. 
Math.  Stat.,  vol.  20,  pp.  477-501,  1949. 

[3]  D.  Chazan,  M.  Zakai,  and  J.  Ziv,  “Improved  lower  bounds  on 
signal  parameter  estimation,”  IEEE  Trans.  Info.  Theory,  vol. 
IT-21,  pp.  90-93,  Jan.  1975. 
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Similarly,  the  transition  from  (25)  to  (26)  starts  at  the  3dB 
lower  point  rj  =  j  where 


■Kt/I'  = 

and  completes  at  77  =  5  where  the  third  term  in  (20)  is  equal  to  the 
first  term,  i.e.. 


i(^  +  ^  = _ - _ 

2V2Sw’  5{ul  +  W^/l2y 


(29) 


Accordingly,  the  lower  bound  for  UWB  signals  has  five  re¬ 
gions,  with  two  thresholds  (the  second  and  fourth  lines)  between 
three  disjointed  segments  (the  first,  third  and  fifth  lines)  for  differ¬ 
ent  77  =  WT  ■  SNR, 


e2  > 


DV12 

threshold 

1(22L  1  -'2 

4  ^  2V2W77  ' 

threshold 

_ TT _ 

r,(wl  +  W^/12) 


rj  <  a 
a  <  rj  <  P 
P  <  rj  <  y 
7  <  77  <  5 
S  <V, 


(30) 


with  threshold  expressed  in  (20)  as  the  sum  of  four  terms. 


4.  SIMULATION 

The  UWB  spectrum  defined  by  the  FCC  is  used  in  a  simulation, 
with  center  frequency  6.85GIIz,  bandwidth  7.5GHz  and  EIRP  emis¬ 
sion  -4L3dBm.  Fig.  2  compares  our  simplified  ZZB,  the  ZZB  ob¬ 
tained  numerically,  and  the  MLE  performance.  (For  this  case,  the 
WWB  is  not  applicable.)  The  simulation  shows  that  our  simplified 
ZZB  is  close  to  both  the  numerically  obtained  ZZB,  and  the  MLE. 
Our  bound  is  divided  into  five  segments,  with  two  thresholds  di¬ 
viding  the  WT -SNR  domain  into  three  disjointed  segments.  When 
WT  ■  SNR  <  a,  the  envelope-related  ambiguities  are  dominant 
and  the  MSE  bound  only  depends  on  the  expected  maximum  time 
delay  D  as  112.  Since  WT l2Tt  1,  when  WT  ■  SNR  <  1, 

that  is,  when  SNR  — >  0,  the  observed  signals  are  completely  dom¬ 
inated  by  noise  and  nearly  useless  for  time  delay  estimation,  i.e., 
the  performance  is  equivalent  to  making  a  random  decision.  When 
/3  <  kUTSNR  <  7,  the  lower  bound  is  contributed  by  both  the  un¬ 
resolved  phase  ambiguity  and  envelope  ambiguity.  /?  is  a  function 
of  both  tJo,  W  and  D.  7  is  defined  as  a  3dB  lower  point  where 
the  MSE  is  half  of  the  one  at  /3.  Obviously,  ft  and  7  increase 
with  D.  The  lower  bound  matches  quite  well  with  the  CRB  when 
WT  ■  SNR  is  very  large  (linear  region  for  WT  ■  SNR  >  20).  This 
region  is  noise  dominated.  Both  phase  and  envelope  ambiguities 
can  be  resolved  completely. 
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Fig.  2.  Comparison  of  the  simplified  bound  with  the  ZZB,  and 
MLE. 
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